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Abstract 



We consider a dilute homogenous atomic Bose-Einstein condensate with two 
non-degenerate internal energy levels. We discuss the case in which the two 
components achieve a state of chemical equilibrium in the presence of an ex- 
ternal radiation field which couples the two states. The presence of the radia- 
tion field can result in new ground states for the condensate as a consequence 
of the lowering of the condensate energy due to the interaction energy with 
the field. We analyze the ground state energy as a function of the coupling 
constants for the two-body interactions, the Rabi frequency of the radiation 
field, and the detuning of the field. We also give explicit expressions for the 
quasiparticle excitation spectrum. 



I. INTRODUCTION 

The recent experimental realization of trapped Bose-Einstein condensates (BEC) with 
internal degrees of freedom [Q] has sparked much theoretical and experimental study of the 
properties of multicomponent condensates. Of fundamental importance is the structure of 
the ground state and the energy spectrum of collective excitations above the condensate 
ground state. This has been explored extensively for the case of spinor condensates in 
which the internal degrees of freedom correspond to the different Zeeman states of a par- 
ticular hyperfine manifold such as the F = 1 manifold in 23 Na. The two-body interaction 
Hamiltonian for a spinor condensate is invariant under rotations in spin space, since for 
two-body collisions, the s-wave scattering length can depend only on the total spin angular 
momentum of the two atoms owing to the rotational symmetry of the collision. In con- 
trast, for two component condensates, such as 87 Rb, in which the two states correspond to 
non-degenerate internal states of the atoms (i.e. states which differ in either their principal 
quantum number or total angular momentum quantum numbers), the two-body interaction 
is not symmetric with respect to pseudo-spin SU(2) transformations since the s-wave scatter- 
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ing length depends on the internal eigenstates of the two atoms. If an external field couples 
the two internal states, one would expect effects that are not present in spinor condensates. 

The ground state and energy spectrum of the quasiparticle excitations of a two compo- 
nent homogenous condensate in the presence of a coupling field has been calculated in Ref. 
||. The coupling field allows atoms to make transitions between the two internal states. 
In that article, as well as in the calculations to be presented below, it is assumed that the 
two components of the condensate are in a state of chemical equilibrium (i. e. the chemical 
potential for the two components are equal 0]). This condition is equivalent to assuming 
that the total number of atoms in the condensate is fixed rather than the number of atoms in 
each of the two components. For nonzero temperature, the relative number of atoms in the 
two components is determined by the condition that the atoms are in thermal equilibrium 

However, for zero temperature (which is the case considered in Ref. [[|), the chemical 
potential may be identified with the energy of a single atom in the condensate. In this 
case, chemical equilibrium corresponds to the condensate being in a stationary state of the 
system. The stationary states of a free atom interacting with an external radiation field 
are known as dressed states in quantum optics ||. Consequently, chemical equilibrium for 
a two-component condensate corresponds to the direct generalization of the atomic dressed 
states f7]. 

Under the condition of chemical equilibrium, Goldstein and Meystre found that the 
quasiparticle spectrum could contain states having imaginary energies, indicating the onset 
of instabilities in the condensate. This is a rather surprising and unexpected result; the 
energy spectrum for quasiparticle excitations above the ground state of an interacting many 
body system are expected to be real. Given these results, we have calculated the ground 
state and quasiparticle spectrum for the model of a two component condensate discussed in 
0. We obtain an expression for the ground state energy density as a function of the relative 
concentration of the two components of the condensate, given a fixed total number of atoms 
in the condensate. For a coupling field of finite strength, the energy density can exhibit 
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both maxima and minima as a function of the number of atoms in one of the components. 
The minima correspond to points of stable equilibrium. As one might expect, the excitation 
spectra calculated around these new minima are real for all momenta. For values of the 
relative concentration other than those corresponding to the minima, in the presence of 
relaxation, the system will be driven to a point of stable equilibrium. The instabilities 
found in Ref. fU can be traced to an expansion about a point of unstable equilibrium. 

In section II, the Hamiltonian for a two component condensate, expressed in a form 
which emphasizes the broken SU(2) symmetry, is used to study the ground state energy 
as a function of the relative fraction of atoms in each of the two states. In section III, the 
quasiparticle spectrum is derived. 



II. HAMILTONIAN AND GROUND STATE ENERGY 



The Hamiltonian operator for a two component BEC may be written as 



H = H 1 + H 2 



(1) 



where Hi is the single body Hamiltonian given by 



Ht = J d 3 r Ufl(r 
Ml 



# (r)+*i(r 



h 2 v 2 
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+ V b (r) 
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#t(r)¥ 6 (r) + #5(r)* a (r) 



(2) 



and H 2 is the two-body interaction, 

H 2 = \j d 3 rd 3 r' {fr a (r)& a (r')U a (r - r')# (r')*a(r) 



^(r)^(rO^(r - rO^(r')^(r) 



(3) 
(4) 



The operators ^(r) are the annihilation operators for an atom in state i = {a, b} at 



position r which satisfy the Bosonic commutation relation, 
with all other commutators being zero. There is a spatially uniform field which couples the 
two states and has a Rabi frequency given by TZ > 0. The atomic field operators have been 



written in a field interaction representation which is rotating at the frequency of the external 
field, LU e . Consequently, there appears in Eq. (g) the detuning, 5 = oj — u e , where hu Q is 
the energy difference between the two internal states. In Eq. (Q), the Vi(r) are external 
potentials which we take to be zero since we wish to consider only the case of a homogenous 
condensate. 

For the two-body interaction, we use a contact potential, Ui(r — r') = Vi5(r — r'). The 
coupling constants, Vi, are expressed in terms of the s-wave scattering lengths by Vi = 47r ^ ai 
where a a and are the scattering lengths for collisions between two atoms in states a and b, 
respectively while a x is the scattering length for collisions between atoms in different internal 
states. In accordance with ||, H 2 can be simplified by assuming that V a ~ V& = V s ^ V x . 

It is possible to rewrite H% in a simplified form that emphasizes the lack of SU{2) 
symmetry. To do this one expresses the field operators as two-component spinors, 



(5) 



The two-body Hamiltonian may then be expressed as 



H 2 = -J rf 3 rrfV^ t (r)^ t (r')V / (r, r')^(r')^(r); (6a) 
Vj(r, r') = \{V S + K.)I(r)I(r')5(r - r') + ^(V s - V x )5(r - v')a z {v)a z {v'); (6b) 

where I(r) is the identity matrix and cr^(r) is the Pauli spin matrix which act on the spinor 
at r. The single body Hamiltonian may also be expressed in the compact form, 

Hl = J ciW(r) |-_ I + y^+--,) *W. (7) 

The previously mentioned lack of 577(2) symmetry in H 2 is now obvious. If we carry out 
an active rotation that diagonalizes Hi, H 2 will not be invariant since the interactions are 
proportional to a z {r)a z {r') instead of cr(r) • cr(r'). 

One can now write a c-number function for the ground state energy. This is accomplished 
by rewriting the field operators as ^(r) = 0j(r) + 54?i(r) where <fii(r) is the condensate wave 
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function which is defined as 4>i(r) = ^,;(r)y = yfple lSi and the expectation value is taken 
with respect to the condensate state |J. Here p(r) = p a (r) + Pbij) is the total density 
which is conserved and S a b = S a — Sb is the relative phase between the two components. By 
making the replacement ^(r) — > 4>i(r) in Eq. (|T|), one obtains an energy functional for the 
condensate. For a homogenous stationary system, Bose-Einstein condensation occurs in the 
state with zero momentum and as such, the 0«(r) are independent of r. The energy density 
for the homogenous condensate with volume V is, 

y = y (Pa - Pb) + hTZ^papb cos(S a b) + ^(V s - V x ) (p a - p h f + ^(V s + V x ) p 2 (8a) 

= y (2 Pa - p) + ^ v /p a (p-p a )cos(s afe ) + hy s - v x ) (2 Pa - pf + hy s + v x ) P 2 (8b) 

Notice that Eq. ( |8a| ) may also be derived from i^H^j = E — i (V s (p 2 + pi) + 2V x p a pb) where 
the expectation value is taken with respect to the wave function 

ic> ^(/¥ e " 4 ° + /S e,s,a »)" io> (9) 

where a\ Q is the creation operator for an atom in state i with zero momentum (see 
next section), N oi = p { V, and N = N oa + N ob . Note that ^ (V s {p 2 a + p 2 b ) + 2V x p a p b ) 
is an intensive quantity whereas the total energy, E a , is extensive and consequently 
Y (V s (p 2 + pi) + 2V x p a pb) is negligible in the thermodynamic limit. Finally, it is clear from 
Eq. (§) that S a b = and n correspond to symmetric and antisymmetric superpositions of 
the internal atomic states. 

In the remainder of the paper, we consider the case in which there is a nonvanishing 
coupling between the two components which allows the atoms to make transitions between 
states a and b. In other words, TZ is never identically zero although it may be infinitesimally 
small, 1Z ~ 0. If TZ = 0, then the number of atoms in states a and b would be separately 
conserved and the condensate energy would be given Eq. (^) for fixed p a and p b - However, 
when Bose-Einstein condensation occurs in the presence of a coupling fied, the condensate 
will form in that state which minimizes Eq. (^) for fixed total number of atoms (p =const.); 
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this is the limit we consider. For non-interacting atoms, this would correspond to the lowest 
energy dressed state. 

Consequently, the ground state of the condensate will be a function of p, 5, TZ, S a b, 
and V s — V x . Note that Steel and Collett have carried out a fully quantum mechanical 
calculation of the ground state for small condensates (a few hundred atoms) under the same 
conditions Similarly, the ground state of a one dimensional inhomogenous condensate 
has been studied under these conditions by Blakie et al. using numerical solutions of the 
time independent coupled Gross-Pitaevskii equations 0. In both these cases, however, the 
relative roles played by the mean-field interactions and interaction with the external field in 
determining the ground state of the condensate is not as physically clear as it is for the case 
of a homogenous condensate. 

In the following two subsections we classify the extrema of Eq. (^a|) for the two cases 
of TZ ~ and 6 = 0. The key point is that nonzero (V s — V x ) can significantly modify the 
ground state structure from what one would expect based on the single body Hamiltonian. 
However, before proceeding it is helpful to make a few definitions. We define the polarization 
of the condensate to be £ = \p a — pb\ / p. Consequently, a polarized condensate corresponds 
to £ = 1 and an unpolarized condensate would correspond to £ = while < £ < 1 
represents a state of partial polarization. 



A. 71 ~ 

In this case we set TZ equal to zero in Eq. (|8a|). When 5 = 0, the ground state will either 
be unpolarized for V s > V x or polarized for V s < V x . When 5 > and V s < V x , the ground 
state corresponding to the minimum of E D /V occurs at pb = p, while E D /V is a maximum 

a t Pa = 2(v - v ) ^ Pl^" ^ or ^ e case > Vx (again with 5 > 0), p a = 2 (v- 5 v ) + P/^ now 
corresponds to the minimum for E Q /V. Notice that when TZ ~ 0, the relative phase, S a b, 
is arbitrary. For 5 < 0, the results are the same if one interchanges states a and b. Note 
that for binary condensates, where N a and N b are seperately conserved, V s < V x would lead 
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to a phase separation of the condensates into two components that occupy nonoverlapping 
regions in space |L3] |T3j fll^f . 



B. 5 = 

When TZ is finite, the ground state can exhibit interesting new structure. It is easy to 
see from Eq. (paf), that the energy is an extremum only if Sab — 0, n. By requiring that 
E a /V be an extremum with respect to p a , one finds from d(E Q /V) /dp a , that the extrema 
are located at p a = p/2 and p a = p^ where 

Pi ±} = ~ (l ± V 1 " V(2p«) 2 ) (10) 

and 

Note that the extrema at pffi are degenerate and occur only when \pa\ > 1/2 and 
a I cos(S' a fe) > 0. There are four cases to consider: 

(i) pa < 1/2 and S ao = 0. In this case the only extremum is at p a = p/2 and this is a 
global maximum of E /V. 

(ii) pa > 1/2 and S a b = 0. In this case p a = p/2 is now a global minimum of E Q /V. The 
extrema at p a = are global maxima of E Q /V. 

(Hi) pa > —1/2 and S ao = vr. Again, p a = p/2 is a global minimum of E /V. 

(iv) pa < —1/2 and S a b = tt- The extremum at p a = p/2 has now become the global 
maximum of E a /V. The extrema at p a = p^ are the global minima in the ground state 
energy. Note that when pa — >• oo, the minima are located at p a = or p a = p. 

The nature of these extrema may be understood by a consideration of the physical 
meaning of the parameter pa. Since the (V s + V x ) term in E Q /V simply gives an overall 
constant energy, it may neglected and, as such, the only relevant mean field interaction 
energy for determining the ground state is the \(V S — V x ) (p a — p ) 2 term. Therefore, pa 



S 



is the ratio of the mean- field energy per atom, ~ (V s — V x )p, to the atom-field interaction 
energy, ~ K1Z. Consequently, for \pa\ > 1/2 the mean-field interactions dominate the ground 
state energy while for \pa\ < 1/2 the energy of the ground state is dominated by the atom- 
field interaction energy. For pa < —1/2 the extremum at p a = p/2 is a maximum, regardless 
of S a b, since the mean-field interactions favor a polarized ground state (£ = 1) in this limit. 
Similarly, for pa > 1/2, p a = p/2 is always a minimum since the mean-field interactions 
favor an unpolarized state (£ = 0). 

On the other hand, for —1/2 < pa < 1/2, the single body atom-field interaction dom- 
inates over the mean-field interactions. In this limit, the minimum of E Q /N = p~ 1 (E /V) 
coincides with the lowest energy dressed state of the free atoms. For 5 = 0, the lowest energy 
dressed is the antisymmetric state ^ (\a) — \b)) with energy —hlZ/2 which corresponds to 
p a = p/2 and S a b = tt. Consequently, the state p a = p/2 and S ao = vr is the global minimum 
of E Q /V in the interval —1/2 < pa < 1/2. The other dressed state for 5 = is the symmet- 
ric state 4| (|a) + with energy hlZ/2 which is not the ground state of the system. The 
symmetric state corresponds to p a = p/2 and S a b = which is a global maximum in the 
interval -1/2 < pa < 1/2. 

For case (ii) and (iv), the energy difference between the extremum at p/2 and p^ is 
given by 



AE/V = i (E(p a = p/2) - E{ Pa = p a +y 



phK 



- pa + cos(Sab) 



;i2) 



4pa 

The energy difference helps to elucidate the transition of the central extremum at \pa\ = 1/2. 
One can see that as \pa\ — > 1/2 from above, p)^> — > p/2 and AE/V — > so that the three 
extrema merge at \pa\ = 1/2 and for \pa\ < 1/2, there is a single extremum at p a = p/2. 
This is illustrated in Figures 1 and 2. 

One should note that E /V always possess at least one global minimum for finite 1Z 
but this minimum does not necessarily correspond to the state with equal population in 
the two components 0]. It is interesting to note that for pa < —1/2 the condensate state 
exhibits another broken symmetry in addition to the usual broken U(l) gauge symmetry 
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since condensation occurs at either or p^ 



III. QUASIPARTICLE SPECTRUM 

In this section the Bogoliubov prescription is used to linearize the Hamiltonian around 
the ground states discussed in the previous section. A canonical transformation is then used 
to diagonalize the Hamiltonian and find the spectrum of elementary excitations above the 
condensate. 

At this point it is advantageous to introduce the grand canonical Hamiltonian, K = 
H - fjiN. Here, N = N a + N b = J d 3 r fi&t (r)* a (r) + ^J(r)^ 6 (r)) is the total number 
operator and p is the chemical potential. We have assumed that the system is in a state of 
chemical equilibrium so that p = p a = p& where p, = J^- is the chemical potential of the 
two components. The chemical potential insures that the ground state expectation value, 
K^j, is a minimum with respect to the total number of atoms. The chemical potential may 
also be interpreted as a Lagrange multiplier which insures that I^N^j is conserved [p}| . A 
pair of equations for p may be derived by requiring that E [4> a , </>&]— p J d 3 r (J0 a | 2 + |0f>| 2 ) 
be an extremum with respect to the variations 5(f)* a and 5<pl, 

^4>a = Y</>a + + (V s \(j)a\ 2 + V x |0fe| 2 ) (j) a ] (13a) 

i4 h = + + (y s \<p b \ 2 + v x |0 a | 2 ) k (i3b) 



It should be noted that Eqs. (|13a[ - |13b|) are equivalent to the condition that (K\ is an 



extremum with respect to the total number of atoms which is found by varying N oa while 
keeping N ob fixed and vice versa. In contrast, the extrema for E D found in section II cor- 
responded to finding (^f 2 ")^ = 0, i.e. the extrema for variations in N oa for a fixed total 



number of atoms in the condensate. As such, Eqs. ( |13a| - |13b| ) serve to define the chemical 



potential, but not the relative fraction of <p a and </>&, which is determined by minimizing Eq. 

For a homogenous condensate we can expand ^(r) in a basis of plane wave states, 
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*i(r) = ^ E P ai P e ip ' r/n , with & = ^ (a i0 ) and 
following expression for K, 



8ij5p,p>- This gives the 



^ = 2^ | + T ~ a " paap + ^2m ~ T ~ M J G ^ a6p + ~2~ V a «p Gb P + 

"'"otT 53 {^s (^api^ap2^ a P3^ap4 + ^bpi^6p 2 ^ fc P3^bp4) + 2V^.a)j pi dJp 2 a a p 3 a5p 4 | ; (14) 

Pl+P2=P3+P4 

which may be linearized around the ground state solutions found in section II by making the 
replacement a i0 — > (a i0 ) and keeping only the lowest order quadratic terms in the operators 
for p 7^ 0. By utilizing Eqs. ( |13a| - |13b| ), the resulting expression for K is 

K = E - pN Q +> h p a Vs — J — COs6 a 6 a' O ap + h PfeKs — J — COsS a 

1 \2m 2 y p a J \2m 2 y p fe 



p 2 M 



Tift 



Q'fepO'ap 



p^O 



^ap^a— p 



+ + V x y/p a p b cosSa^j (a[ p a bp + aj p a ap ) + -p^ (a^ p al_ p + 
+ ^P&Vs (4 P 4-p + «bp«6- P ) + V x y/p a p b cosS ab (al p al_ p + a ap a&- P 
A Hamiltonian which is quadratic in bosonic operators with the general form, 



1 - ' < ',t At 1 



H = H + - J2 Aj(p)aJ p a}_ p + - £ ^(p^p^-p + £ ^(P^ 



(16) 



«,j,P^O iJiP^O i,j,p¥=o 

where -Bjj(p) is Hermitian and Ajj(p) is a symmetric matrix may be diagonalized by a 
canonical transformation. This is done by defining quasiparticle annihilation and creation 
operators, i>i p and b\ , respectively, given by 

(17) 



i 

which satisfy bosonic commutation relations. The quasiparticle operators obey the equation 
of motion 



k P ,H 



(18) 



which is consistent with H = J2i, P ^o £ i(pWi p hp + E vac . Using Eqs. ( [18] ) and ( |T?D it is easy 



(15) 



to show that the quasiparticle energies satisfy the eigenvalue equation [ |I~2|j , 
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£j(p)uk(p) = E ( B ij(p) u ij(p) + A ij{p)vij(p)) ; 

j 

-£,(pW(p) = E (^(P)%(P) + ^,(pMi(p)) ! (19a) 

j 

Since £/(p) are not the eigenvalues of a Hermitian matrix, there is no guarantee that they 
will be real. However, when the £/(p) are complex, the i>i p and b\ do not satisfy bosonic 
commutation relations and consequently, the quasiparticles may no longer be interpreted as 
bosons. In addition, if (u^(p), %(p)) are a solution with eigenvalue £/(p), then (t> ; * (p), u^(p)) 
are a solution with eigenvalue — £i(p). However, only Ei(p) > are physically significant since 
the energy of the system must be bounded from below. 

Since Eq. ( jT5|) has the form of Eq. (|T6|), one may directly apply Eqs. ( |19a| ) to calculate 



the energy spectrum of the quasiparticle excitations above the condensate. As in the previous 
section, we focus on the two cases of 1Z ~ and 5 = and limit the discussion to the case 
V x >V S (a < 0). 

A. K ~ 

For 5 > 0, the ground state is given by pb = p and p a = 0. In this case Eq. ( |15|) has the 
form 



k = e - pN + 2 (f^U»p + £ b(p)S p 4J 



(20) 



where £_e(p) = y^(^ + 2pV^^ and 5 P is the quasiparticle annihilation operator for 
excitations in state b and is related to db p through the canonical transformation, db p = 
cosh<y2 p i?p — sinh^p^lp [|ll|]. For 5 < 0, condensation occurs in p a = p and Eq. (p0|) re- 
mains valid with the interchange of a ap <-> a^p. Finally, for 5 = and V x > V s , condensation 
only occurs in the state pb = p or p a = p since the two states are degenerate. Bose conden- 
sation will not occur in the state pb = p a = p/2 as given in Ref. || since this represents a 
maximum in E Q /V. 
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B. 5 = 0, K > 



For pa < and 5 = 0, the only minima in E Q /V occur for S a b = tt. For the extremum at 
p a = p/2, the quasiparticle energies are 



MPJ 



MP) 



2m 



Ml+p(V s -V x ) 



(21a) 
(21b) 



which agrees with the results obtained in Ref. ||. It is possible for e+(p) to be imaginary 
and this situation corresponds to the instabilities mentioned in [[J. The condition that £+(p) 
be real for all momenta is pa > —1/2 which is exactly the condition that p a = p/2 be a 
minimum of E D /V. Consequently, Eqs. ( 21a - 21b| ) can be interpreted as the quasiparticle 
spectrum for excitations above the condensate only when pa > —1/2. When Eqs. ( 21a - [21b ) 
are extended to pa < —1/2, they no longer correspond to physically meaningful results 
since e±(p) then correspond to expansions of small amplitude oscillations around an energy 
maximum. 

When pa < —1/2, the location of the minima of E Q /V are p a = p^>. Since Eq. ( |T5D is 
symmetric under the interchange of a and b atoms, the quasiparticle spectrum will be the 
same at p a + ^ and p a ~^- The excitation spectrum about p a ~^ (or p a + ^) is significantly more 
complicated than at p a = p/2 and is given by 



—1 



p" / pr 



2 mK 2 m +pV ) + \^ 2 ^- 1 
(Ml) 4 ((2pa) 2 - l) 2 + AihTZf ((2pa) 2 - l) 

1 \ , 



(22a) 
(22b) 

m (2 P a) +P V,) + L + ^\^ 



\2m 



(pVsY + 1 + 



2pa / 



2pa)(KR)(KR-2pV s 



(22c) 



Note that the ± in e± (p) correspond to the two branches of the excitation spectrum for each 
of p a + ^ and p a ~^ and not to p a + ^ and p a ~^ individually. It is obvious that Ei > and one can 
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prove that S2 > for 2pa < — 1. One may also show that Sf — S2/4 > which implies that 
£±(p) will be real for 2pa < — 1. Thus the quasiparticle spectrum is real as expected. 
When p = 0, e_(0) = and e+{0) = 



p 2 (V x — V s ) — (T1R) 2 which shows that one 
branch is gapless while the other branch contains a gap for pa < —1/2. It is easy to show 
that when pa = —1/2, Eqs. ( |22a| - |22c| ) agree with Eqs. ( |21a| - |21bD which indicates that the 
quasiparticle spectrum varies continuously with pa. There are several important limiting 
cases for Eq. ( |22a| ). First one may consider the limit that pa — > —1/2 while ^-,pV s ~ hTZ. 
In this case the excitation spectrum reduces to 



Mp) 



\ 2m (^2m ^ ^ x s ' 



(23) 



Note that Eq. fl23|) is not valid in the long wavelength phonon regime since j- ~ H1Z ~ 
p {Vx — Vs) ■ Another simple case is when the two minima approach the edges at p a = 

2 

and p a = p which corresponds to the condition \pa\ 3> 1/2 while at the same time #- <^ 
H1Z, pV x ,pV s (long wavelength limit). In this case the two branches of Eq. ( 22a ) are given 
by 



MP 



MP) 



IpVs 



m 



-p; 



m 



P 



p{Vx ~ Vs) + 2m~- 



(24a) 

(24b) 
(24c) 



Notice that in this case, the e- (p) branch corresponds to phonons with a speed of sound given 
by u = \J ^ that is independent of the interspecies scattering and that £+(p) corresponds 
to single particle excitations with a mean-field shift in the energy. 

2 

Finally, one may consider the limit that pa — > —1/2 and hIZ which corresponds 

to the long wavelength limit in which the two minima are displaced only slightly from p/2. 
In this case the excitation spectrum has the form 



e±(p) 



Si ± 



{KKY {{2paf - 1) - 2pV,£) + ^{WZf ({2paf - 1) (25) 
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which for pV s ^> H1Z, e±(p) simplifies even further to 

f v I p(Vx + Vs) (oa , 

£ ' (P) = V 2m P] (26a) 

(26b) 

Again, e_ (p) corresponds to phonon excitations with a speed of sound which depends on the 
average of V s and V x . However, the e+(p) branch can not be given a simple interpretation 
in terms of either phonon like collective excitations or single particle excitations since both 
the p 2 term and the term in brackets are of comparable magnitude. 

IV. DISCUSSION 

In this paper we have analyzed the ground state energy of a homogeonous two-component 
Bose-Einstein condensate interacting with a spatially uniform radiation field which couples 
the two internal states of the condensates. We have argued that Bose-Einstein condensation 
occurs in the state which minimizes E Q as a function of the density of one the components 
(since the total density is fixed) and the relative phase between the two components. When 
the proper ground state for the condensate is chosen, the quasiparticle excitations above the 
ground state have real energies and no instabilities occur. The fact that the quasiparticle 
spectrum is real and exhibit no instabilities may not come as a surprise, but these results 
differ from those obtained in Ref. 0. The instabilities found in Ref. || originate from an 
expansion of small amplitude collective excitations about a point of unstable equilibrium, 
represented by the global maximum of Eq. fl8af) . 

Throughout this paper it has been assumed that the condensate is in a state of chemical 
equilibrium (i.e. \i a = fa = P where fa is the chemical potential for component i). As such 
the ground state energy is found by minimizing E a subject to the constraint p a + Pb = P- An 
interesting extension of this work would be to consider the case when the two components are 
in a state of thermal equilibrium such that Bose condensation can occur, but not in a state 
of chemical equilibrium such as in [[J. In this case the population of the two components 
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is fixed by some external pumping and decay mechanism and the ground state energy of 
the condensate will be given by Eq. (|8a| ) for fixed p a and p^. In this case one might expect 
instabilities to occur in the excitation spectrum since the condensate state is not in general 
the minimum of Eq. (|8a|). Normally, these instabilities would indicate a transition to a 



phase separated state PT5[] . However, if the pumping and decay are spatially uniform and 
time independent so that the relative densities at each point in space are fixed externally, 
then it will be impossible for phase separation to occur. The dynamics of the condensate in 
this case would be interesting to explore. 
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Figure 1. Plot of E Q /V for Sat, = and pa > showing the transition from a central 
maximum at p a = p/2 for pa < 1/2 to a minimum for pa > 1/2 along with the 
development of maxima at p a = p^ + \ p^~\ Note that E Q /V is in units of phlZ and we 
have subtracted off the constant term \(V S + V x )p 2 . 

Figure 2. Plot of E Q /V for S a f, = ir and pa < showing the transition from a central 
minimum at p a = p/2 for pa > —1/2 to a maximum for pa < —1/2 along with the 
development of minima at p a = p^ + \ p(~\ Note that E Q /V is in units of pMZ and we 
have subtracted off the constant term \{V S + V x )p 2 
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